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INVISCID HYPERSONIC PLOW OVER A BLUNT BODY WITH 


HIGH RATES OF MASS AND HEAT TRANSFER 
By E. Dale Martin 
Ames Research Center 


SUMMARY 


A study of steady hypersonic flow over a "blunt Body with a very high rate 
of continuous mass transfer from the body surface, with heat conduction, is 
presented. An approximate analysis is made of the inviscid flow over a sphere 
with large mass flux. Lighthill's well-known constant-density inviscid blunt- 
body flow solution was previously extended by Vinokur and Sanders and by 
Cresci and Libby to include inviscid flow at a high rate out of the body sur- 
face. The present analysis further extends Cresci and Libby 1 s solution to 
include the effects of the boundary shock wave, a thin viscous region domi- 
nated by viscous-compressive stresses (rather than viscous shear) and the 
associated heat conduction in the gas adjacent to the body surface. These 
effects are important to consider when there is a significant source of trans- 
lational nonequilibrium at the body surface (such as absorption of intense 
radiation and the accompanying large heat conduction) . 


INTRODUCTION 


Hypersonic flow over a blunt body with a high rate of mass transfer has 
been studied by a number of investigators. For example, in an early investi- 
gation, Hoshizaki (ref. l) made numerical calculations from the incompressible 
Navier-Stokes equations to determine the effects of mass transfer on the heat- 
transfer rate to blunt bodies. Vinokur and Sanders (ref. 2 ) extended the 
constant-density inviscid solution of Lighthill (ref. 3 ) to include mass 
injection from the body. The calculations of reference 2 were made for the 
very special case wherein the injection is vectored so that the flow from the 
body is irrotational . Libby (ref. 4), Cresci and Libby (ref. 5 ) , and Fox and 
Libby (ref. 6) also studied blunt-body flows with large mass-injection rates. 
Of particular interest here is the appendix of reference 5^ which contains a 
solution to the inviscid, incompressible-flow equations that is equivalent 
(except for the vectored injection) to that given in reference 2 . The invis- 
cid solutions of Vinokur and Sanders and of Cresci and Libby are further 
discussed in references 7 and 8. 

The magnitude of the mass flux from the body determines the nature of the 
interaction with the outer flow and the nature of the effects on the body 
motion and heating. When the mass flux is small enough, the outer flow is 
not affected appreciably, and boundary^layer theory is applicable if the 
Reynolds number p^V^do/ |i 2 is large. (Notation is defined in appendix A.) 

It will be shown in a later section that the bulk of the flow in the inner 
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region is essentially inviscid if the mass flux (or, more appropriately, Re^) 
is large; that a thick viscous region exists only for moderate values of mass 
flux; and that a boundary layer exists only for small mass flux (small Re^). 

A thirj. region of flow dominated by viscous-compressive stresses (Trr) and 
heat conduction is expected to occur, in general, at a boundary surface from 
•which a gas flows at a rapid rate with significant heat conduction at the sur- 
face (e.g., because of absorption of intense radiation there), as discussed 
in references 9 and 10. Thus, even though the bulk of the flow from the body 
is inviscid, a thin layer immediately adjacent to the body generally can have 
significant viscous effects which are quite different from the viscous- 
shearing effects (e.g., significant T rcp) that characterize a boundary layer. 
In this region, designated as a boundary shock wave, the density, temperature, 
pressure, and velocity may vary rapidly. In the inviscid limit, the boundary 
shock wave is of vanishing thickness, and the flow variables change discon- 
tinuously across it (in the same sense as across a shock wave). Significant 
heat transfer at the wall and large Re^ are necessary and sufficient condi- 
tions for the presence of the boundary shock wave (ref. 10) whether the 
efflux is supersonic or subsonic. If the heat transfer is small, the effects 
of the boundary shock wave are negligible. If the efflux is supersonic rela- 
tive to the surface (Mb > l), the limiting case of the boundary shock wave as 
the heat conduction approaches zero is a simple detached shock wave near the 
surface . 

The effects of the boundary shock wave did not appear in the analyses of 
references 1, 2, 4, and 6. Although the viscous terms were included in 
the incompressible form of the Navier-Stokes equations used in reference 1, 
the boundary- shock effects would not appear in a purely constant -density 
analysis unless discontinuities were allowed in order to account for rapid 
variations in density, normal velocity, etc. Similarly, the inviscid-flow 
analyses of reference 2 and the appendix of reference 5 contained neither 
viscous-compressive effects (effects of T rr ) nor discontinuous flow variables 
to account for possible viscous effects at the body. Those analyses also 
would not apply to supersonic efflux because, with no heat conduction at the 
wall, there would have to be a detached shock wave somewhere between the body 
surface and the stagnation point . Boundary- shock effects did not appear in 
the boundary-layer analyses of references k and 6 because the terms in the 
Navier-Stokes equations that can be important very near the surface when the 
mass flow normal to the surface is large (terms containing T rr) are omitted 
in the boundary-layer equations used in those analyses. Since the boundary- 
shock effects are not important in some cases, the results of the previous 
investigations are valid and accurate (under appropriate conditions) even 
though account was not taken of a possible boundary shock wave. 

It is hoped that the present study will shed some light on the nature of 
the division of the flow into inviscid and thin viscous regions. After justi- 
fication of the approach used, the flow is analyzed by treating constant- 
density inviscid regions with lines of discontinuity to represent thin viscous 
regions. (Treating the internal structure of the viscous regions is beyond 
the scope of this report, but rigorous derivation of the boundary conditions 
for the inviscid regions, based on recognition of the presence of thin viscous 
regions, is contained in an appendix.) The analysis includes the effects of 
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the boundary shock wave, which makes the flow calculation correct in the 
extension to flows with more severe conditions than can be treated correctly 
by the previous analyses. The inviscid-flow solution is the same as that in 
the appendix of reference 5 except in the application of boundary conditions 
at the body surface. That is, the inviscid solution is the same as Cresci 
and Libby* s if their boundary parameters are the same as values outside a 
boundary shock in the present analysis. Thus, one main purpose of this report 
is to relate the inviscid solution to appropriate boundary parameters under 
the conditions where the boundary shock wave is expected to occur according to 
the theory of references 9 and 10 . The approach presented here may therefore 
be considered to be an extension of the analyses in reference 2 and the 
appendix of reference 5 to cases with significant heat conduction at the 
boundary surface due to a source of translational nonequilibrium, such as 
intense radiation being absorbed there. 


DIVISION OF FLOW INTO INVISCID AND THIN VISCOUS REGIONS 


The flow of air and of gas emanating from the body is assumed to be 
governed by the Navier-Stokes equations. They may be expressed symbolically 
as 


convective terms = — • (viscous terms) 

Re 


(i) 


The convective terms and the viscous terms are made dimensionless with respect 
to the appropriate flow quantities in each region of interest so that the 
Reynolds number for the general case is 


Re 


p ref V ref l ref 
^ref 


( 2 ) 


where the reference quantities in a given region are constant values 
representing actual significant flow quantities in that region. 

It is well known that, if the Reynolds number is large, the right side of 
equation (l) is negligible, and the flow is said to be inviscid, except in 
very thin regions where the viscous terms, which contain higher order deriva- 
tives of velocities than the convective terms, are large. In the limit as 
Re -> 00 ^ it is only in regions for which the thickness is of the order of some 
power of (Re) -1 that these viscous terms are important. In the inviscid 
limit the thin regions are represented as discontinuities in the flow vari- 
ables that are changing rapidly there (the discontinuities being either at a 
boundary or at an internal location, where either a boundary condition or 
internal continuity of the solution cannot be satisfied because the order of 
the differential equation, or system of equations, is reduced). The structure 
of the flow inside these ”quick-transition regions” (cf. refs. 11 - 13 ) could 
be analyzed by making boundary-layer-type transformations on the equations and 
making a uniformly valid approximate calculation by dropping out the terms of 
higher order in the small parameter (Re)"" 1 . 
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Consider first the case of no mass transfer (V-b = 0) and refer to 
figure l(a) , which represents schematically the flow in the forward region of 
a blunt body that is axially symmetric with respect to the free-stream direc- 
tion. A significant parameter determining both the thickness of the shock 
wave, 6 S w; and. the boundary-layer thickness, &bZ> in. comparison to the 
thickness of the shock layer, d D , is 


Re 2 


V-z 


(3) 


where subscript 2 indicates conditions on the axis immediately behind the 
shock wave and where p, p, and V are, respectively, the density, viscosity, 
and velocity there. As is well known. 


e S 



(v 2 - Voo j 0)1 


b bZ 


% l A /_ - 1 / 2 . 
= = 0 (Re„ J 



l as Re 2 00 


w 


in general. Thus, when = 0 and Re 2 is large, the flow of air in front 
of the body (fig. 1(a)) is inviscid except for the thin quick-transition 
regions of the shock wave and the boundary layer (shaded areas). 


Now consider the cases where mass transfer from the body exists (V^ ^ 0), 
and in particular the regimes of flow for which Re 2 » 1 and for which 


Reb 



(5) 


is very small, of order unity, or very large (see fig. l(b), (c), (d)). The 
case of small mass transfer (f^V^ small) is denoted by Re-^ « 1. The 
"intermediate regime" is Re^ = 0(l); and the case of large mass transfer is 
denoted by Re^ » 1. 

From the above discussion we see that Re^ is the determining factor as 
to whether or not the bulk of the flow from the body (i.e., the inner gas 
layer) is free of viscous effects. (Note appendix B.) For very small mass 
transfer (Re^ « l) the inner gas layer is entirely viscous and is simply 
part of the boundary layer (fig. 1(b)). The mass transfer does not change the 
order of magnitude of the boundary-layer thickness if Re^ < 0(\/Re 2 ) . As the 
mass flux increases, say to where Re^ = 0(l), the inner layer becomes a 
thick viscous layer (fig. 1(c)). Mathematically speaking, it is not a bound- 
ary layer because it is not vanishingly thin as Re 2 °o and the requirements 
for the validity of boundary-layer theory are violated, although the viscous 
part of the air in the outer layer does become vanishingly thin as Re 2 -> °o # 

As the mass flux increases still further, so that Re-^ » 1, the inner gas 
layer must become essentially inviscid (fig. 1(d)), leaving viscous effects 
confined to the interface layer and the boundary shock wave at the body 
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surface. The interface layer is essentially a "blown-off” boundary layer 
that has carried with it the viscous-shearing effects (effects of T rC p) that 
characterized the boundary layer (fig* 1(b)). The viscous-compressive effects 
(effects of T rr ) that became of the same order as the viscous shearing 
effects in the thick viscous layer (fig* 1 (c)) are left confined to the 
boundary shock wave (fig* 1(d)). Note that no part of the flow in the bound- 
ary shock wave need be supersonic if there is heat transfer in the gas at the 
boundary (ref. 9 ). Furthermore, the boundary shock wave behaves in the same 
manner with varying Re^ as does a shock wave with varying Re 2 (see ref. 9 ). 
Its thickness for large Re^ is of order ?h/Pb v b^ that 


^ _ ^bsw 

€ bsw - — ^ — 
d i 


OCRe^ 1 ) 


as Re^-, -»• 00 
with V e - V b / 0 


(6) 


where the viscosity coefficient p is of the same order as the shear- 
viscosity coefficient p. The viscosity coefficients are related by 


M- - 



K 


( 7 ) 


where k is the bulk viscosity coefficient (see ref. l4, p. 337 )* One would 
expect the order of magnitude of the relative thickness of the interface 
layer , , to be given by 

= hi = 0(Re 2 /2 ) + 0 (Re^ 1/2 ) as Re 2 00 and Re^ -*■ °o ( 8 ) 

if it can be treated like a boundary layer . However , this should probably not 
be assumed a priori without a careful analysis. It would be true if the tan- 
gential velocity at the interface is discontinuous in the inviscid solution. 


Interface Boundary 

Shock wave layer shock wave 



Sketch (a) 


With the assumption that the 
entire outer flow is governed by 
Re 2 and the inner flow by Reh^ 
the limiting case for which both 
Re 2 and Ret, -> <» will leave two 
inviscid regions with possible dis- 
continuities at their boundaries: 
the shock wave , the interface , and 
the body. Whether or not discon- 
tinuities will exist there in the 
inviscid limit will depend on satis- 
faction of the conservation equa- 
tions across these vanishingly thin 
regions, to be considered below. 

Sketch (a) may aid the under- 
standing of flow containing a bound- 
ary shock wave . Thermodynamic 
pressure, p, and mass density, p. 
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are plotted qualitatively across the several regions of flow on the axial 
streamline when Re 2 and Re^ are hoth large and when some source of signifi- 
cant translational nonequilibrium is present at r^ • Note that in measure- 
ments of pressure, there would he no way of physically separating the 
thermodynamic pressure, p, from the normal viscous stress in the radial direc- 
tion, For small enough pV-^ 2 , the total normal stress in the radial 

direction, p + (-T rr ), does not vary signif icantly across the boundary shock, 
although the density, p, does if a strong boundary shock is present (i.e., 
with significant T rr ) . 

ASSUMPTIONS AND APPROXIMATIONS FOR FLOW CALCULATION 


For the flow calculation it is assumed, first of all, that Re 2 and Re-^ 
are large enough that both the air flow in the outer layer and the flow 
originating from the body are essentially inviscid. It is assumed therefore 
that all viscous effects are confined to very thin layers at the shock wave, 
the interface, and at the body surface, as described above. These very thin 
quick-transition regions are represented by surfaces of discontinuity in the 
inviscid analysis. The discontinuous quantities are those that vary rapidly 
inside the thin viscous regions, except that the density at the interface is 
generally discontinuous also in the viscous flow, if diffusion is neglected. 

It is assumed that either the interface layer is stable near the axis of 
symmetry or the effects of instability are negligible (see appendix B ) . 

Only the sphere will be considered for the body shape . In order to 
obtain a separable solution to the equations, the shock wave and the inter- 
face will be assumed to be concentric with the body, and the normal velocity 
of mass flow from the body (just outside the boundary shock wave) will vary 
as the cosine of the body angle cp. These assumptions are automatically 
included if we simply assume: (l) the body has constant curvature, and ( 2 ) 

the flow solution is separable. That is, the latter assumption is equivalent 
to specifying all the boundary conditions in such a way as to be compatible 
with a separable solution to the equations. 

For simplicity, the density is taken to be constant in each of the invis- 
cid regions, but it may vary across the surfaces of discontinuity represent- 
ing rapid transitions. (The density varies rapidly across the boundary shock, 
as well as across the shock wave.) The outer inviscid solution will then be 
identical to Lighthill's constant-density solution for the sphere (ref. 3 ) 
(with the interface as the equivalent body in that solution), since the con- 
ditions across the shock wave are the same . The solution for the flow between 
the interface and the body will be different from those of Vinokur and 
Sanders (ref. 2) and Cresci and Libby (ref. appendix), for the same bound- 
ary parameters, because of the presence of the boundary shock wave. (For the 
same values of parameters outside the boundary shock wave, the inviscid solu- 
tion is the same as given in the appendix of reference 5 • The relationship 
of the inviscid solution to the boundary parameters, under conditions where 
the boundary-shock-wave theory applies (ref. 10 ), is of most concern here.) 

The assumptions of constant density outside the boundary shock and the sepa- 
rable flow solution simply require that the mass flux from the body, as well 
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as the normal velocity outside the boundary shock, he specified as varying as 
the cosine of the body angle <p. These conditions will be evident when the 
solution is formulated. 

Purely radial velocity of the mass flow out of the body surface will be 
assumed. The flow through the boundary shock wave will then be normal to the 
surface. In the solution of Vinokur and Sanders (ref. 2) the mass injection 
was assumed to be vectored to make the flow irrotational . The vorticity, or 
rotationality, of the flow from the body depends simply on the distribution 
of the components of velocity at the surface, which can be specified indepen- 
dently of the flow problem, since they depend entirely on the injection (or 
ablation) process . 

The limiting case of very small heat conduction (let) 0) in a super- 
sonic efflux of gas at the boundary, for which the boundary shock becomes a 
detached shock wave (ref. 9)* will not be treated here. Thus we assume 
<lc-b / 0 lf Mb > 1- 

The mass flux from the body may be either subsonic or supersonic for a 
boundary shock wave to occur (see refs. 9 a *icL 10) . We do not restrict M 5 
except to require / 0 , since = 0 would be generally incompatible with 
the limit Re-^ 00 • 

^ A value of unity is assumed for the longitudinal Prandtl number, 

Pr = pcp/k c , that is, the Prandtl number based on the viscosity coefficient 
[x that occurs in the purely viscous-compressive stress, T rr = p.(dv/dr) , in 
the boundary shock wave. With that assumption the simple relationships 
developed in reference 9 lor th e conditions across the boundary shock wave 
can be used. Liepmann and his co-workers have pointed out (ref. 15) that 
Pr « 1 for most real gases. 


INVISCID CONSERVATION EQUATIONS AND THE BOUNDARY CONDITIONS 
FOR THE INVISCID REGIONS 


Inviscid Conservation Equations 

For steady inviscid flow the equations of conservation of mass, conserva- 
tion of momentum (Euler equations), and conservation of energy are: 


div(pV) = 0 


(9) 


p(V • grad)V + grad p = 0 


(10) 


(ID 

with which appropriate equations of state and a relation for q r must be 
used. We take for the equations of state in the inviscid regions, as 
described above. 


p(V • grad)^h + — ■ V 2 ^ + div q r = 0 
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p = p = constants (r^ < r < r^ + dj_) 

p = p + = constant, (iv + dj_ < r < + d) 



Sketch (b) 


Equations ( 9 ), (10), and (ll) maybe 
written conveniently in axisymmetric 
spherical coordinates (r,cp) for the 
sphere (see sketch (b)). 

A transformation that simplifies 
both the equations and their solution 
(but involves no approximation) is 
the following: 


£ = log(r/r b ) 


for which 




(13a) 


(13b) 


In terms of dimensionless quantities 
u, v, p, p, and k, equation (9) and 
the components of ( 10 ) become 


(pue 2 ^ sin cp) + (pve 2 ^ sin cp) ^ = 0 

pCuu^ + vUj. + uv) = -kp c 
p(uv m + w f - u 2 ) = -kp 


( 15 a) 

(151) 


With dimensionless enthalpy h, equation (ll) becomes 


pufh + ~ u 2 + v 2 J + pv^h + ^ u 2 + v 2 J 


1 „ 2 \ “ kr div 


p V 3 

Hx> V 00 


In equations (l4) to (l 6 ), the subscripts cp and £ denote partial differen- 
tiation. (in these equations p is still variable.) The equations of state 
for the inviscid regions become 


p = p 


p /p + = cr 2 = constant for 0 < £ < 


P = P = 1 


for 
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where 


= log ( r i/r^, ) = log(l + di/i^) = log D ± 
C s = iogCrg/rt,) = log(l + d/i^) = log D 


(18) 


It is convenient to introduce a dimensionless Stokes stream function \|/, 
defined by 


. - 2t . 

i = pue b sin cp ; 


— 2. t 

^cp = -nve b 


pve b sin cp 


( 19 ) 


which automatically satisfies equation (l 4 ) . Equations ( 15 ) may be combined 
by using equations (17) and by cross-differentiating and equating to 

to obtain a nonlinear partial differential equation for which applies 
for 0 = £ e < £ < ^ or for < £ < £ s . Because constant density is 
assumed for the inviscid portions of the flow, the energy equation is 
uncoupled from the equations for conservation of mass and momentum and, hence, 
need not be considered for the inviscid-flow calculation (except in the con- 
ditions across the discontinuities to determine p” and p + ) . 


Boundary Conditions 

Since the flows on both sides of the shock wave and on both sides of the 
interface layer are inviscid, the inviscid-flow equations for the conservation 
principles can be used to determine conditions across (cf. appendix C) . How- 
ever, the boundary shock wave, like a boundary layer, has inviscid flow only 
on one side, so the conservation equations of viscous flow must be used to 
find conditions across the boundary shock wave. This has been done (refs. 9 
and 10 ) for Pr = 1 and for a plane boundary surface. For a curved boundary 
surface and a very thin boundary shock wave (i.e., to first order) with flow 
normal to it, the results will be the same. However, the flow quantities will 
vary along the boundary shock wave . 

The conditions behind the shock are (cf. appendix C): 
a "k £ = £s : 


-k cos cp 

(20a) 

sin cp 

(20b) 

(l - k)cos 2 cp 

(20c) 

i (l - k 2 )cos 2 cp 

(20d) 


where k = p^/ p 2 is determined by an equation of state (or equilibrium-air 
tables) in combination with equations ( 20 ) (using cp = 0 ) . It will be 
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convenient to use (20c) in the following form (making use also of eq. ( 15 a)): 

£ = , -kp^ = uu cp + vu ^ + uv = +2k(l - k)cos cp sin cp (2l) 

The interface (£ = is defined as the surface between the region of 
gas flow that has come through the shock wave and the region of gas flow 
emanating from the body. This definition can be represented by 

£ = i ± > * = 0 (22a) 

which, with the concentric interface (i.e., £± = constant, not a function of 
9), also implies 


t = £i , v = 0 


(22b) 


The conditions to be applied at the inviscid interface are (cf. appendix C): 
at £ = £ ± ■ 


v = 0 


P = P + 


P ( uu <p) = p (uu ffi ) 


<P ' 


(23a) 

(23b) 

(23c) 


At the body (in the inviscid flow, outside the boundary shock wave) it 
will be found necessary to specify: p = p e , v=v e , u = Ug (where 

p e = p = P - /p + = a 2 = constant) in order that the inviscid solution be deter- 
mined. (These quantities are equivalent to those "at the body" in the appen- 
dix of ref. 5 .) When a boundary shock is expected, these quantities must be 
related to conditions at the boundary, within the viscous boundary-shock-wave 
region. The viscous-flow equations must be used to determine the conditions 
across the boundary shock, whose relative thickness (ref. 9) is 0(Re^ x ) . 

For the calculation of results, one would like to specify values at the 
boundary, say M^, T^, and the mass flux PfeV^, then determine the conditions 
outside the boundary shock (p e , v e , Ug), which would be used as boundary con- 
ditions for the inviscid solution. However, it will be found most convenient 
for the calculation to specify values of a parameter N containing P e V eo 2 . 

As will be seen, then, p eo will be determined from the solution to the equa- 
tions so that Mg 0 2 = Pe^eo 2 /^Peo w iH be known. If = Mj-, ( cp) and T^q are 
also specified, then the complete inviscid solution, the conditions across the 
boundary shock, and values of quantities of interest at the boundary, includ- 
ing the mass flux at cp = 0, (PbVt>) 0 , w iH be determined. Thus, specifying 
M)-,, T-b 0 , and P e V eo 2 will give solutions corresponding to Mb, Tb 0 , and some 
value of (PbVb) 0 * lb bas been shown (refs. 9 and 10 ) that, for Pr = 1, per- 
fect gas in the boundary shock, and q r = constant across, and for known 
values of Mb and M e , the conditions across the boundary shock wave for flow 
normal to the layer (in present nomenclature) are: 
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v e /vb = a 

p e / Pb = 1 / a 

T e /T b = 1 + (1/2) (7 - 1)^(1 - a 2 ) = a 2 ly^ 2 /M g 


( 24 ) 


where 


a = 


Me 

Mb 


1 + (1/2) (7 - l)Mb 
.1 + (1/2) (7 - l)Me 2 J 


2 -1 1/2 


( 25 ) 


and the heat -conduct ion parameter for conduction in the gas at the boundary is 


c h, * 


-Icb 


(1/2) PbV 


= 2 r/i + -i_v 6- + ^V 

L \ 7M e 2 / \ 7Mb /- 


(26) 


To these conditions we must add for normal injection , 1 

u e = ■% = 0 


( 27 ) 


In general, the quantities found from the above conditions are variable over 
cp. The respective variations will be specified as required to produce a 
separable solution for the inviscid flow, as discussed above. For a separable 
solution it is necessary that v e vary as cos cp and at most as sin cp: 


at £ = C = 0 : 


v = v e = v eo cos cp = 


V, 


eo 


V„ 


cos cp 


u = u e = u eo sin cp = 0 


( 28 ) 


Once the ratios in equations (24) are known for all cp (with specifica- 
tion of M^, TbQ, and p e V eQ 2 and with p e found from the inviscid solution), 
the following quantities can be calculated in order: 


1 ln a private communication. Prof. N • Rott (Univ. of Calif, at Los 
Angeles) has pointed out that if a significant tangential velocity component 

is present in the boundary shock wave, its treatment would be analogous to 
the "shock slip" condition for the tangential velocity component in an oblique 
shock wave (e.g., see ref. 1 6 ). This condition would have to be considered if 
the mass injection is "vectored" as suggested in reference 2 . 
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p e (<p) 

P e /Pb 


P b (9) = 




Tb(<p) 




(29) 


2 

Different values of P e V eo could be specified initially and the solution 
iterated until the desired mass flux, p^V-^ , is obtained. 

It should be noted that when the mass transfer is a result of ablation 
of the body surface caused by high radiative heating from the hot layer of 
air between the shock wave and the interface, the specification of arbitrary 
values of M^, T b , and p^V-^ could be replaced by three appropriate equations 
similar to those used in references 9 and 10 for the one-dimensional 
case: a vaporization-rate equation for , Clapeyron’s equation (or "modi- 

fied Clapeyron's eq.") for T]-, = T^p^,) , and an equation relating Pb^ to 
the radiative heat flux. Then the radiation is assumed to be absorbed within 
an infinitely thin layer just inside the body surface; some of that heat is 
conducted into the body, so the heat conduction into the body from the surface 
is not the same as the conduction in the vapor just outside the surface. 


SOLUTION OF THE EQUATIONS FOR THE INVISCID REGIONS 


Separation of Variables 

The partial differential equation for f obtained from equations (15), 
(17), and (19) with the boundary conditions (20), (21), (23), and (28) has 
the separable solution 
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i|f = sin 2 cp f(£) 

where the differential equation for f ( £ ) is 


(30) 


f.r. _ 5 f" + 2f • + 8 f = 0 , (0 = 5 e < C < 5 ± or q < S < S s ) 

where ( )' = d( )/d£. The boundary conditions become: 
at 5 = £ s = log D : 

f = kD 2 /^ 

f = d 2 r 

f"' = (D 2 /k) (1 - k + 2k 2 ) j 


(31) 


(32) 


at S = £. + = log D ± : 
at (= q~ = log D ± : 


f = 0 


(33) 


f = 0 

(f) 2 = [(f’) 2 3‘ = (^) [(f') 2 ] + 

\P / 

at g = £ e = 0 : 



— P Ueo ® 


(3*) 


(35) 


Equation (3l) is a third-order linear ordinary differential equation (with 
constant coefficients because of the transformation ( 13 ); its solution will 
have three unknown constants of integration for the region 5i < £ < £ s and 
three unknown constants of integration for the region 0 < £ < In addi- 

tion, and £ (or and D) are unknown, making a total of eight unknown 
constants to be evaluated by the eight boundary conditions in equations ( 32 ) 
through (35) • 


Equations for Pressure, Velocity, and Vorticity 

The boundary conditions not yet used, (20c) and ( 23 b), will determine the 
constant of integration for the pressure in each of the two regions . When 
equations ( 19 ) and ( 30 ) are used, equations ( 15 ) become 
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-kpp 9 = (sin cp cos cp)e 4 ^[(f') 2 - 2ff" + 2ff»] 


(36a) 

(36b) 

(37) 


-kpp^ = (sin 2 cp)e 4 5[2ff* - (f 1 ) 2 ] + (cos 2 cp)e 4 5[4ff' - 8s 2 ] 
which may be integrated to obtain 

-kpp = (1/2) (sin 2 cp)e _4 k(f') 2 - 2ff" + 2ff'] + e _4 ^(2f 2 ) + C 
where the conditions (20c) and (23b) now determine that 

c = c + = | k 2 - k for < 5 < £s 

= C“ = p"C + = a 2 (i k 2 - k) for 0 < £ < 

Outside the boundary shock wave, conditions (35) in (37) give (for Uqq = 0): 


(38) 


P. = 1 - ^ k - 


1 ,, (^eo) 2 ^eo f "(°) 


2k 


2k 


sin^ cp 


(39) 


In terms of the function f(£), the dimensionless velocity components 


are: 


u = 


sin cp f '(0 

- 2 £ 
pe b 


(ho) 


V = 


-2 cos cp f(£) 

pe^ 


and the dimensionless vorticity is 

— > 

r-i curl V 


w = 


= 


e ^(Vcp - U - u^) 


(41) 


sin cp 

pe 3 5 


(2f + f ' - f") 
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Integration of the Differential Equation 
and Calculation of Results 


The solution to equation ( 3 l) nia y be written 

f(5) = f + (0 = c£e 4( ^i } + c +e 2(S " 5l) + cje~ ( ^ i} 
= e*(r/rp 4 + c^(r/r ± ) 2 + c^r/rU' 1 

(Ci < £ < 5 S ) 


(tea) 


and 


fU) = f-(5) = c-e 4( ^~ Cl) + c^e 2(C ^ + c^e (5 ^ 


4 b __ 

= c~(r/r ± ) + c"(r/ri ) 2 + c 3 (r/r 1 ) 


-1 


( 42 b) 


(0 < 5 < q) J 

The constants c l7 c+, c 4 , c~, c“, c~, D, and are evaluated by substi- 
tuting conditions (32) through (35) into the solution (te). The results are 


C+ _ (1 - k) 2 

4 = 4 k - 1 . 

C3 (l - k) (l - 6k)D 3 

D 2 10 kD 2 

D 2 6k 

D 2 15k 


where 

- = D_ = r b + d = *s 
~ D i ^ + d ± "" r ± 

is the appropriate root of 

3(1 - k) 2 + 5(4k - 1)D 2 + 2(1 - k)(l - 6k) D 5 = 0 

Also obtained are: 

^ -6>,= + 2D,* - 

> 


Per 


/ 2B 


^2 

Pa 


-D 4 

1 


k + 


cm. 


■eo 


D- 2 ) — 
0 1 J 2B 


- 2 - = ( - 2 D.- + 4+2 

Pa 


CTU eo r\ 2 \ D i 


D ± j — 
1 / 2B 


(* 3 ) 

( 44 ) 


(te) 


( 46 ) 
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where Dj_ is the appropriate root of 


and where 


root of 


A(Dq) 

N - 0 

(47) 

B(D ± ) 

CVJ 

O C\J 
<L> 

> > 

a 4 ’ a 

II 

o 

6*“ 

III 

(48) 

V 00 oo 


s p-/p + = p" = P e = kp e / Poo 

(*9) 


P = P(k) = 


_ u(£ = Zp _ = H) 


sin cp 


sin cp 


j/c+X /cj\ _ /c+\ 

~ \pf/ + %y ' VDfy 




15k 


■ 6( - + 5(4k - 1) - D 3 (l - k)(l - 6k) 


(50) 


and 


A(Dj.) = [3 - 2(aueo/P)]Di 5 - 5D ± 3 - (dUgo/P^ 2 + 2 + 3((n^ 0 /p) (51) 

B(Di) = Dj 5 + 5D ± 2 - 6 (52) 

For small dj/r^ (D^ l), equation (47) may he expanded (for = 0) as: 

r h 3 W 3 W 


+ • • • 


(53) 


and the inverse expansion, conveniently used for evaluating di/r^ for given 
small N, is 


— = N - | N 2 + ii N 3 - +• 
3 9 


(ueo = 0) 


(54) 


It is particularly useful to note that f + ( t,) /Dj_ 2 depends only on k 
(since c 1 + /D± 2 , c 2 + /D^ 2 , and c 3 + / D i 2 are functions of k and D, and 
5 = D(k)), so that (l/Dj. 2 )f + and its derivatives can he calculated versus £ 
independently of the solution for 0 < £ < £i* Similarly, f - (£)/Pcf depends 
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only on the parameter N = ov eo /p for u eo = 0 (since c^/Pa, c 2 “/Pcr, and 
c 3 - /p a are functions only of Di for \Xq Q = 0, and Dj_ = Dj_(N) for u^ Q = 0 ), 
and can he computed independently of the solution for < £ < £ s . 

Another result of special interest here is that, for a given flight con- 
dition (given Voo, p^, p^) , for which k and P are known, f" (£ = 0 )/a is 
known so that p e (eq. ( 39 )) depends only on crv eo or on N (eq. (48)); 
hence, p e can he found for given N regardless of the value of either a or 
v eo • Then the conditions across the boundary shock wave can he worked out in 
the manner outlined in the section "Boundary Conditions," hy specifying values 
of M^, Tb 0 , and P e V eo 2 or N . (See eqs . (24) through ( 29 ).) Equivalently 
(and often more conveniently) one could specify different values of D-^ 

(rather than N), for which N is easily calculated from equation (47). 


DISCUSSION OF RESULTS 


The ratio of the width of the outer inviscid region to the radius of the 
interface, shown in figure 2(a), is identical to that given hy Vinokur and 
Sanders (ref. 2) and hy Cresci and Lihhy (ref. 5 ); values of that ratio are 
also identical to the values given hy Lighthill's solution (ref. 3) for the 
ratio of the shock standoff distance to the hody radius for the same values 
of k = P 0 q /p 2 in the case of no mass injection. Shown in figure 2(h) is the 
ratio of the thickness of the inner inviscid region to the hody radius versus 
the parameter N = ov eo /p. Note that this "inviscid blowing parameter" N 
is based on conditions outside the boundary shock wave, and not on conditions 
immediately adjacent to the hody surface when a boundary shock is present. 

The dimensionless tangential velocity on either side of the interface 
(Uj_ + or Uj_“) can he found from the function P plotted versus k = Poo/ P 2 in 
figure 3. This function P is identical to the dimensionless stagnation- 
point velocity gradient from Lighthill's solution. 

Figure 4 shows examples of the tangential -velocity profile for the 
particular cases for which k = 0 .06 and dj_/rb =0.3. Note that the func- 
tion plotted in figure 4(a) is continuous at the interface, hut that the 
velocity itself is discontinuous if p~/p + is not unity (see fig. 4(h), in 
which p“/p + =2). 

The remaining calculations, for quantities at or near the hody surface 
(figs. 5 to ll), were all made for = 20. 83 km/sec, altitude = 60 km (for 
which k = 0.06554 and P = 0.4190). The ratio of specific heats in the gas 
emanating from the hody was taken to he 7 = 7/5- To illustrate the variation 
of quantities along the boundary shock wave (i.e., with varying cp) in fig- 
ures 5> 6, and 7j the particular value 0.025 was arbitrarily chosen for the 
parameter N. Figure 5 shows the variation of p e , P e V e 2 , and M e . Note that 
these variations correspond to the assumed cosine distribution for the veloc- 
ity outside the boundary shock wave. For figures 6 and 7 it was also neces- 
sary to choose a value for Mb- For these calculations, Mb = l/3 was 
arbitrarily chosen. The values for Tb 0 an( i the molecular weight m need 
not he specified for these results. 
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For the same flight condition (V m = 20 .83 km/sec, altitude = 60 km), the 
values of various quantities of interest relating to the boundary shock wave 
in the gas on the axis at or very near the body surface are shown plotted 
versus the inviscid blowing parameter N (for 7 = 7/5 in the inner region) in 
figures 8 through 11. One need not specify either Mb, Tbo> or 111 in fig- 
ure 8 , but calculating the ratios across the boundary shock, conditions at 
the body surface, and certain conditions outside the boundary shock wave 
requires that Mb be specified (figs. 9 to ll) . It is pertinent to note 
again (as in the discussion of sketch (a) ) regarding the pressure variation 
across the boundary shock, that a measured pressure normal to the r direc- 
tion would be the total normal stress, p + (-T rr ), rather than the purely 
thermodynamic pressure p . Note also that the heat conduction parameter Cj lc 
(fig. 10 (a)) expresses the heat conduction in the gas at the boundary; its 
relationship to the heat conduction inside the body surface depends on the 
mechanism of the mass transfer (refer to refs. 9 and 10) . For all the 
quantities plotted in these figures, neither Tbo nor m need be specified. 
Worthy of particular notice is the result shown in figure 10 (d), the relation- 
ship between the mass flux (p^Vb) and the parameter N referred to in above 
discussion. 


If one makes the approximation 



(55) 


one can plot the results given in figures 5 through 11 in more general form. 
Then, the solution for appropriate dimensionless quantities would depend only 
on k, rather than on the entire flight condition with specified altitude and 
flight velocity. For example, with the flight condition used in figure for 
which P V ^ 2 = 1.533X10 6 dynes/cm 2 . 


and 


PeV 


10 3 dynes/cm 2 


= 1533 


p eV 

p V 2 
■ 00 00 


10 s dynes/cm 2 


= 1-533 


Pe 


PnX 


(56) 


so that, with the approximation ( 55 ), the values on the three curves on fig- 
ure 5 are equivalent, respectively, to the two quantities on the right sides 
of equations ( 56 ) and Mg for just k = 0.06554, N = 0.025, and 7 = 7 / 5 . 
Similarly (with use of the approximation (55) and with the flight condition 
(altitude and velocity) replaced by k = 0.06554), the results of figure 6 
apply directly, and those on figure 7 for Cj^/lOO and T^/T^p apply directly. 


whereas the values on the curve labeled 


50 gm/mole 
m 



become 


18 



the values of 


/50 gm/mole \ /Tb 0 \ /l. 533X10 6 cm 2 /sec 2 \ /Pb\ 

V m ) \1°V \ V 2 J VV 

N ✓ \ V x * CO ' X ' 


The remaining 


results of figures 7 to 11 can he worked out for the more general condition 
in a similar manner (with, e.g., the quantities Me, p e /p b , P e /%> T e/ T b> 

C hc applying directly) • Of special interest is a further result, plotted in 
figure 12: the dimensionless heat-conduction flux in the gas at the body 

versus another "blowing parameter," based on conditions at the boundary. 


F = 


/Pbo^bo 2 \ 

V 4oV/ 


(57) 


showing values on various lines of constant N or lines of constant M 5 • For 
given N and M^, the blowing parameter is calculated from 


p V ^ 
_bbo 


p V * 

*oo 00 




( 58 ) 


and the dimensionless heat conduction is found in the form 





( 59 ) 


Obviously, the variation of the heat transfer with the blowing parameter 
depends on the simultaneous variation of Mb* For example, if Mb remained 

constant as F increased (where Mb = Vp^V^ 2 (f/ 7Pbo ) ) j then 

( V^/RTbo ) 17 2 1 -q.c b / ( 1/2 ) P^Voo 3 ] would decrease. As pointed out above (para- 
graph following eqs . ( 29 )), the complete determination depends on other con- 
ditions of the problem, including the mechanism of the mass transfer . 

In discussing the above results it must be recognized that the boundary- 
shock-wave theory, and the inviscid-flow solution, are only valid asymptoti- 
cally as Re b 00 - Hence, for sufficiently small N or F, the results given 
here are not valid. As N and F approach zero, the inner flow becomes fully 
viscous and then approaches a boundary-layer-type flow, so that, in the limit, 
the ratio P e /Pb must go to unity; the heat conduction becomes the usual 
convective heat transfer of boundary-layer theory; the normal viscous stress 
becomes insignificant in comparison to the thermodynamic pressure p; and the 
viscous shear, T rC p, becomes significant. (However, even a value of N equal 
to 0.02 or 0.01 can correspond to a substantially high value of Re^ • ) As an 
example, if for sufficiently large N a constant value of Mb were appro- 
priate as N varies (e.g., Mb ** l/3> as estimated in refs. 9 and 10), then 
as N 0, the variation of P e /Pb would be expected to follow a curve such 
as the qualitative dashed curve sketched on figure 9(a) - Similarly, the value 
of C bc in figure 10(a) would first approach a curve of variation appropriate 
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to "moderately large" values of mass flux, and eventually to that appropriate 
to boundary-layer theory as N and Re-^ approach zero. 


Ames Research Center 

National Aeronautics and Space Administration 

Moffett Field, California 9^°35> Aug. 9, 1967 
12^07-02-23-00-21 
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APPENDIX A 


C 

Ch c 

c p> c v 


Ci>C 2 ,c 3 


D 



F 

f 

h 

h 

k 

k c 

M 


PRINCIPAL NOTATION 

constant of integration in pressure; evaluated in equation ( 38 ) 

heat-conduction parameter, equation ( 26 ) 

specific heats at constant pressure and constant volume, 
respectively 

constants of integration; evaluated in equations (43) and (46) 

1 + A.fa 

r b 


(eq. (44)) 


1 + — = — 


d o + d i 

thickness of the layer of gas emanating from the body 

thickness of the layer of air between the shock wave and the 
interface 

Pu V. Z\l/2. 
hdo bo 

p V 2 

nx> 00 


function of £ in the stream function, equation ( 30 ) 
specific enthalpy 
h - h™ 


V 2 

CO 


shock density ratio, — 


coefficient of thermal conductivity 

V 


Mach number. 


speed of sound 
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m 


N 


Pi- 


molecular weight 

ov /> P ^ 2x ;L/,2 

u eo _ / k *e eo 


3 VP 2 p V 2 

X *00 00 

pLCp ^ 

(Prandtl number based on p) 
*c 


P 

P 

q. 

— > 

q. 

R 

Re l 

Re- 


thermodynamic pressure 

p - p 


p V 
H oo 00 


magnitude of heat flux 
heat flux vector 


gas constant 


Pb V h d : 


’(?) 


times universal gas constant 


— ) , Reynolds number 




p V d Q 

— — — , Reynolds number 
distance from the sphere center 


u 


temperature 

dimensionless tangential velocity, X tangential velocity 

V 00 


V magnitude of velocity 

V velocity vector 

v dimensionless radial velocity, X radial velocity 

V oo 

Z magnified independent variable in thin region; product of £ 

reciprocal of some small dimension, depending on the region; 
appendix C 

a — (eqs. (2k) and ( 25 )) 


and 

see 
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0 


7 

&bZ 

6 bsw 

5 iZ 

S sw 

«bZ 

€ bsw 


i 

K 

P 

P 

P 

P 


a 

T rcp> T rr 


9 


co 


function of k (eqs. ( 50 )); equivalent to stagnation-point 
velocity gradient in Lighthill solution 

C P 

ratio of specific heats , — , in the gas from the body 

c v 

boundary-layer thickness 

thickness of the boundary shock wave 

thickness of the interface layer 

thickness of the detached shock wave in the air 

&b z 

— , equations (4) 
a o 

, equation ( 6 ) 
d i 

5 

— = — , equat ion ( 8 ) 
d 

— ^ - y equations (4) 
d o 

dimensionless independent variable, equation ( 13 a) 
bulk viscosity coefficient 
shear viscosity coefficient 

+ (J) K (eq. (7)) 

mass density 
P 



respectively viscous-shear stress and viscous-compressive stress 
in the radial direction 

angle between a radial line through a point in the flow and the 
body-free-stream axis 

dimensionless Stokes stream function, defined by equations ( 19 ) 
dimensionless vorticity, equation (4l) 
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Subscripts 


b value just outside the body surface 
c heat conduction 

e value just outside the boundary shock wave 
eo value outside the boundary shock at cp = 0 
i value at the interface 

o value at cp - 0 

r radiative heat transfer (except in T rC p and T rr ) 

s value behind the shock wave 

£ partial derivative with respect to £ 

cp partial derivative with respect to cp (except in T rC p) 

a value behind the shock wave at cp = 0 

00 value in the free stream 

Superscripts 

+ value in the outer air layer, between the interface and the shock wave 

value in the inner gas layer, between the boundary shock wave and the 
interface 


2k 
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APPENDIX B 

ON THE POSSIBILITY OF STANDING EDDIES OR AN UNSTABLE 

INTERFACE LAYER 


We may consider the possibility of regions internal to the flow whose 
essential character is not governed by either Re 2 or Re^, . (See appendix A 
for notation.) The steady flow immediately behind the shock wave is certainly 
governed by Re 2 and that adjacent to the body is certainly governed by Reb • 
But suppose there is some factor at the interface not accounted for that would 
allow one or more regions (standing eddies) of dimension, say, Z*, to develop 
in which typical values of p, V, and p are p*, V*, and p*. Then a Reynolds 
number 

Be* - <2^ (Bl) 

could actually be the governing parameter in that region; and if, for example, 
Re* were of order unity, then there would be a region of viscous flow not of 
vanishing thickness as Re^ and Re 2 oo. Such an event is considered to be 
unlikely, and is simply assumed not to occur when the flow calculation is 
made in the text of the report . 

There is another related possibility that can be eliminated only after 
careful consideration: instability of the interface and the formation of 

eddies . The interface between the air and the gas flowing from the body is a 
contact surface of discontinuity. In the inviscid limit, in the general case, 
the tangential velocity is discontinuous (shown in appendix C) and the inter- 
face layer is a "vortex sheet." According to Rouse (ref. 17, p. 303 ); "Since 
the time of Helmholtz it has been -understood that a vortex sheet is inherently 
unstable and will degenerate into a series of vorticity concentrations if the 
slightest disturbance is present." Prandtl (ref. l 8 , pp. 50-53) discussed sev- 
eral examples of surfaces of discontinuity, all of which break up into eddies. 
In Prandtl f s words: "Owing to fluctuations in the flow the surface of separa- 

tion may take on a slightly wavy form .... The waves advance with a veloc- 
ity which is equal to the mean of the two velocities . . . ." Such a surface 
of discontinuity then "... breaks down into separate eddies" (ref. l 8 , 
pp* 50, 5l)* Prandtl described further the mechanism of the instability. 
According to Prandtl’ s descriptions one might expect the interface (for both 
Re 2 and Re^ very large) to look like the purely hypothetical case depicted 
in sketch (c). Perhaps some physical reasoning can be added that will help 
to explain the development into eddies in the other examples and will also 
bring out a contrasting situation in the present problem. In all the examples 
discussed by Prandtl and in all others (at least in subsonic flow) known to 
the present writer, the pressure in the region where a contact surface is 
breaking up Into eddies is always increasing with streamwise distance. When 
a contact surface of discontinuity is formed in the flow over a sharp edge , 
the breakup (instability) appears to begin at the point where one would expect 
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the pressure to begin increasing 
(where the streamlines begin diverg- 
ing 1 ) . Hence the contact surface is 
unstable in an "adverse pressure 
gradient." The contact surface of 
discontinuity is really a very thin 
viscous layer, like a boundary 
layer, between two flows at high 
Reynolds number. A boundary layer 
tends to separate from its boundary 
in the presence of an adverse pres- 
sure gradient (see discussions by 
Schlichting, ref. 19, pp. 23 and 
26-32, and Prandtl, ref. l 8 , p. I37), 
which is always associated with the 
formation of vortices (ref. 19; 
p. 23). Since the "internal bound- 
ary layer" at the surface of "dis- 
continuity" does not have a wall to 
separate from in the presence of the 
adverse pressure gradient, the vorticity concentrations tend to separate from 
each other and so to form eddies in the same way as a separating boundary 
layer does. Thus it is conjectured that the same phenomenon causes the insta- 
bility and collapse of a contact surface that causes separation of a boundary 
layer, which only occurs in the presence of an adverse pressure gradient. It 
is further conjectured, then, that the Helmholtz instability and eddy forma- 
tion may not occur at a contact surface in the presence of a favorable pres- 
sure gradient. In the present problem, the streamlines are known to be 
converging and the pressure decreasing. Therefore the eddy formation and 
collapse of the contact surface are not expected to occur in this problem, 
especially near the axis of symmetry. 

Since it has not been proved here that the instability will not occur, 
some further remarks may be made. First, it is known that the tendency toward 
instability and eddy formation is greater the higher the Reynolds number. 
Hence, the Reynolds number Re-^ may be large enough that asymptotic solutions 
for high Re-b are useful, but Re-k may not be too large to prevent viscous 
damping of disturbances in the flow pattern that would initiate the eddy for- 
mation (cf. ref. 17 , p. 303 ). That is, the viscous interface layer may be 
thick enough that an unstable "discontinuity configuration" of the flow is not 
too closely approached. 

Second, even if Re^ is very large and if the instability and eddy- 
formation condition can develop in spite of the above conjectures, the region 
in which the eddies occur may be thin enough not to influence the inviscid 
flow significantly, but merely to influence the structure of the viscous 
interface layer, at least in the region of interest near the axis of symmetry. 


Shock wave 



Sketch (c).- Hypothetical unstable 
interface . 


1 From conservation of mass, one finds diverging streamlines to have 
decreasing velocity and, hence, from Bernoulli's equation, increasing 
pressure . 
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Third, if the fluid densities in the two invise id regions are the same, 
then there is no discontinuity in the inviscid velocities and hence no 
likelihood of an unstable interface . 
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APPENDIX C 


DERIVATION OF INVISCID SHOCK-WAVE AND INTERFACE 
BOUNDARY CONDITIONS 


The main purpose of this appendix is to show a consistent and rigorous 
derivation of the conditions applying across the interface layer. Because of 
the convenience of the direct analogy to a similar derivation of conditions 
across the shock wave, the formal derivation of conditions across the shock 
is shown first. 


To derive the conditions across the shock wave and the interface, con- 
sider the shock to have relative thickness e s and the interface layer to 
have relative thickness €j_ = + e . (See appendix A for notation.) 

These thin quick-transition regions are parallel to constant £ lines; that 
is, the flow quantities vary rapidly over £. For the shock wave define 
Z = £/e s so that AZ = 0(l) across the shock wave. Assume q r is constant 
across the shock wave (i.e., no significant radiation is absorbed or emitted 
within the 'shock) . Then the Navier-Stokes equations for the flow through the 
shock are (cf. eqs . (l4), (15), and (l6)): 

(pv) z + 0(e s ) = ol 


pvu^ + 0(e s ) + viscous terms = 0 
pw z + kp z + 0(e s ) + viscous terms = 0 




(Cl) 





+ viscous terms = OJ 


from which we see that, in the inviscid limit (e s 0; viscous terms -> 0 on 
*both sides). 


pv = constant 
u = constant 
(pv)v + kp = constant 


across the shock 





constant J 


(C2) 


that is, the inviscid conservation equations are satisfied across the dis- 
continuity representing the quick transition region only if equations (C2) are 
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satisfied. These well-known conditions across a shock wave give the 
conditions behind the shock as equations (20) in the text of the report. 


The interface (£ = £j_) is defined as the surface between the region of 
gas flow that has come through the shock wave and the region of gas flow 
emanating from the body. (See eqs. (22) in the text.) The conditions across 
the interface can be derived in a manner analogous to that used for the shock 
wave . Let e-j_ + be the relative thickness of the air part of the interface 
layer (£ > £p and €■?“ be the relative thickness of the part in £ < £^. 

Now let Z = (£ - £i)/ e i + for £ > £ ± and Z = (£ - £ 1 )/e i - for £ < £ ± , 
so that AZ = 0(l) across the interface layer. It can be shown that, in the 
analogy with the boundary layer, if u varies rapidly across the interface 
layer, then v = 0(e-j_) = e-jW inside the viscous layer, where e 

e± = €j_ + = 0 (Re a ) -1 / 2 for £ > £ 3 . and where = 0 (Re - b ) -1 ^ 2 for 

£ < Assuming q r is constant across the interface layer, we can 

represent the Navier-Stokes equations (cf. again eqs. (l4) , (15), and (l 6 )) 
as: 


(pu sin cp)^ + (pW sin cp) z + 
ptuujp + Wu z ) + kp^ + 0(e ± ~) + 0(e ± + ) 
kp z + 0(ei“) + 0(e 1 + ) 


P 






+ 0(s ± ) + 0(e 1 + ) 


0(ei“) + 0(€.j_ + ) = 0 
+ viscous terms = 0 
+ viscous terms = 0 

+ viscous terms = 0 


(C3a) 

(C3*) 

(C3c) 

(C3cL) 


In the inviscid limit ( 0, 0, viscous terms ^ 0 on both sides of 

the interface), equations (C3a) and (C3d) give no information about conditions 
across the discontinuity. Equation (C3a) is automatically satisfied, since 
mass does not cross the interface. Equation (C3d) simply says that 
h + (l/2 )u 2 is constant along the inviscid streamlines immediately adjacent 
to both sides of the interface. Equations (C3c) and (C3*b) show, respectively, 
that the inviscid momentum equations are satisfied across the discontinuity 
representing the interface layer only if 


p = constant across the inviscid interface discontinuity 
p(uU(p + vu^) + kp^ = 0 on each side of the inviscid interface discontinuity 

(C^a) 

Since p is continuous across, from the first equation, and v = 0 on both 
sides (eq. (22b)), the second equation requires 

puu^ = continuous across the inviscid interface (C^b) 

The conditions to be applied at the inviscid interface are therefore 
equations ( 23 ) in the text of the report. 
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Figure 1.- Mass- transfer flow regimes. 


33 







k = 6.06 (P = O.itOa), dj/r^ = 0.3 (PN = crv eo = O.I395), arbitrary 

a 3 = p-/p+ 


Figure 4 .- Examples of tangent ial-velocity profile for sphere. 
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Figure 5 •- Distribution along the boundary shock wave of flow parameters 
outside the boundary shock on a sphere at 60 km altitude; V m = 20. 83 km/i 
(k = 0.0655*0 for N = 0.025, 7=7/5* 




Figure 6.- Distributions along tbe boundary shock wave of conditions across 
the boundary shock on a sphere at 60 km altitude; V ro = 20 .83 km/sec 
(k = 0.0655^) for N = 0.025 and for constant Mu = 1/3 and 7 = 7/5* 
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Figure 7 • - Distributions along the boundary shock wave of conditions in the 
gas at the surface of a sphere at 60 km altitude; = 20. 83 km/sec 
(k = 0 .0655*0 for N = 0 . 025 , = 1 / 3 , 7 = 7/5, and arbitrary Tk 0 and m 
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(a) Pressure ratio. 

Figure 9.- Ratios of flow conditions across the boundary shock wave (at the 
axis) on a sphere at 60 km altitude; V = 20 . 83 km/sec (k = 0.0655^*) for 
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(c) Temperature ratio. 
Figure 9*~ Concluded. 
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(c) Density (for given T^ 0 and m) 
Figure 10 — Continued. 
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(a) Density (for given T^q and m) . 

Figure 11— Values (on the axis) of flow parameters outside 
wave on a sphere at 60 km altitude j = 20.83 km/sec (k 
y = 7 / 5 } various Mb, and arbitrary Tbo and m. 
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Figure 12.- Heat conduction versus blowing parameter for k = 0.0655^; 
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